Compressibility Instability of Interacting Electrons in Bilayer Graphene 
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Using the self-consistent Hartree-Fock approximation, we study the compressibility instability of 
the interacting electrons in bilayer graphene. The chemical potential and the compressibility of the 
electrons can be significantly altered by an energy gap (tunable by external gate voltages) between 
the valence and conduction bands. For zero gap case, we show that the homogeneous system is 
stable. When the gap is finite, the compressibility of the electron system becomes negative at low 
carrier doping concentrations and low temperature. The phase diagram distinguishing the stable and 
unstable regions of a typically gapped system in terms of temperature and doping is also presented. 
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Bilayer graphene has attracted considerable atten- 
tion because of its promised application in electronic 
devices^ - — In contrast to the Dirac fermions in mono- 
layer graphene, the energy bands of the free electrons 
in bilayer graphene are hyperbolic and gapless between 
the valence and conduction bands. Most importantly, 
an energy gap opening between the valence and conduc- 
tion bands can be generated and controlled by external 
gate voltages. At low carrier doping, because the energy- 
momentum dispersion relation around the Fermi level 
is relatively fiat, the Coulomb effect is expected to be 
significant in the interacting electron system of bilayer 
graphene. The Coulomb effect has been studied by a 
number of theoretical works using the Hartree-Fock (HF) 
and random-phase approximations*^— One of the ther- 
modynamic quantities directly reflecting the Coulomb ef- 
fect is the electronic compressibility. Recently, this quan- 
tity has been measured by experiments^— For timely 
coordinating with the experimental measurements, it is 
necessary to theoretically study the combined effect of 
Coulomb interactions and the gap opening in the com- 
pressibility. 

In this work, using the self-consistent HF (SCHF) 
approximation^ we investigate the energy bands, the 
chemical potential and its derivative with respect to the 
carrier density for the interacting electrons in bilayer 
graphene. With these results, the compressibility insta- 
bilities for gapped and ungapped electron systems are ex- 
amined. We show that the compressibility is always pos- 
itive for zero gap case, while for the gapped system the 
compressibility becomes negative at low doping and low 
temperature that implies the system would become inho- 
mogeneous. The phase diagram distinguishing the stable 
and unstable regions in the temperature-density plane is 
presented for a typically gapped homogeneous electron 
system in bilayer graphene. SCHF is a conserving ap- 
proximation satisfying microscopic conservation laws^* 
The conserving conditions are crucial in the many parti- 
cle systems. 

The atomic structure of bilayer graphene is shown in 
Fig. 1. The two sublattices in each layer are denoted 
by A (white) and B (black) atoms, respectively. The 




FIG. 1: Structure of Bernal stacking bilayer graphene. The 
energies of intra- and interlayer electron hopping between the 
A (white) and B (black) atoms are given by t and ti, respec- 
tively. 



interlayer distance is c = 3.34 A 1.4a where a ~ 2.4 
A is the lattice constant of monolayer graphene. The 
energy of electron hopping between the nearest-neighbor 
(nn) carbon atoms in each layer is t ps 2.82 eV^i while 
the interlayer nn hopping is t\ « 0.39 eVi^ 2 . 

The Hamiltonian describing the electrons is given by 

P ij 

where t/^ = (c^, c^ M , c£ j0)2 , c^ 2 ) with p = (k,v,a) 
standing for the momentum k, the valley v and spin 
<7, and Cp a e creating an electron in p state at a lat- 
tice of the I layer (with the top layer labeled as 1), 
rii is the electron density at the site i, and Vij is the 
Coulomb interaction between electrons at sites i and j. 
Here, for the low energy electrons in bilayer graphene 
with interlayer hopping, the matrix h is an extension 5,8 
of the Dirac fermions in monolayer graphene . 23 i 24 By 
denoting k x + ik y = fcexp(i0), h can be written as 
h° = T(cj>)h{k)T\<l>) with T{<j>) = Diag(l,e^, e^, e l2 ^) 
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FIG. 2: (color online) Diagrammatic expressions for (a) the 
self-energy E(fc), (b) the irreducible density-density response 
function, and (c) the vertex correction. The solid line with 
an arrow is the Green's function and the wavy line represents 
the Coulomb interaction. 
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in units of Eq = = a = 1. Here A is the energy 

gap parameter, which is a consequence of the potential 
difference between the top and back gates (attached to 
the top and lower layers, respectively). The form of h° 
reminds us to take the transform tp pa = T((f>)^ pa and 
then work with the scalar momentum k. We here take 
the cutoff of k as k c = 1. Under the SCHF approxima- 
tion to the Coulomb interactions, one obtains an effective 
Hamiltonian 
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with h e ff(k) = h(k) + S(fc) and 



^X;tV(l*-#l)c°sM*t*<W.(4) 



V 

/ a a 2a\ 

a a 

a a 

\2a a a 



(5) 



where v^ v {\k — k'\) is the Coulomb interaction, a is the 
angle between k and k', and "ly^ [with p' = (k',v,a)] 
is the /ith component of Q p > . For electrons in the same 
layer, v^ v {q) = 2ne 2 /eq with e = 4. For interlayer elec- 
tron interactions, v^ LV (q) — 27re 2 exp(— cq)/eq. The self- 
energy £(fc) is diagrammatically shown in Fig. 2 (a) 
in terms of the Green's function. In the numerical cal- 
culations, special care needs to be paid to a logarith- 
mic singularity stemming from the azimuthally integral 
of v^dk - k'\) cos 6^,,. 

H c ff can be diagonalized by using the transformation 
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FIG. 3: (color online) Energy band structure of free electrons 
(dashed lines) and the interacting electrons under the self- 
consistent Hartree-Fock approximation at T = (solid lines) 
for (a) A/So = and 5 = and (b) A/E = 0.02 and 
S= 10" 5 . 



where U{k) is an unitary matrix and V pv is the eigen op- 
erator of a quasiparticle with energy E v {k). The chemical 
potential /i is determined by the doped electron density 
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where the factor 4 comes from the spin and valley degen- 
eracy, and f(E) = l/{exp[(E - n)/T] + 1} is the Fermi 
distribution function. The diagonalization of H e s and 
determination of and p, are carried out by iteration 
until the self-consistency is achieved. 

In Fig. 3, we show the energy band structure of the 
electrons at temperature T = for (a) A = and elec- 
tron doping concentration 8 = \/3a 2 n/8 (doped electrons 
per carbon atom) = and (b) A/Eq — 0.02 (the typi- 
cal magnitude in experiments^) and S = 10~ 5 . The 
solid lines are the energies of the electrons renormalized 
by the Coulomb interactions. The four energy bands of 
free electrons denoted by the d ashed lines are given by 
E°(k) = ±[A 2 + tf/2 + k 2 ± ^t\/A+(i{ + 4A 2 )fc 2 ] 1 /2. 
The upper and lower bands of the free electrons are sym- 
metric about zero energy. For the interacting electrons, 
the energies shift downwards because of the self-energy 
contributions, and the upper and lower bands are not 
symmetric anymore. For the zero-gap case, the zero gap 
is not changed by Coulomb interactions. Notice that 
the contribution from the two lower bands to the self- 
energy is included here. If it is a constant for any doping 
and temperature, it can be subtracted from the begin- 
ning. However, in the SCHF, the valence band itself and 
thereby its contribution to the self-energy vary with dop- 
ing and temperature. Especially, in the zero gap case, 
this variation is significant. The contribution from all 
the two lower bands cannot be considered as a constant 
and subtracted from the beginning. 

Figure 4 shows the chemical potential /i as a function 
of 8 for the system with A/Eq = 0.02 at several temper- 
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FIG. 4: (color online) Chemical potential fi as function of 
doping concentration S for A/Eo — at T/To — 1 (black 
dashed line) and A/E = 0.02 at T/T = (blue dotted 
line), T/T = 0.5 (green dash-dot line) and T/T = 1 (red 
solid line). 



FIG. 5: (color online) (d/j./dn)T as function of doping con- 
centration 5 at T — (blue dash-dot line), T/Tq = 0.5 (green 
dashed line), and T/T = 1 (red solid line) for A = 0. The 
black dotted line is the result of the free electrons at T = 0. 



atures and that of the zero-gap system at T/Tq — 1 with 
To = 0.01-Eo ~ 283 K (close to the room temperature). 
For the gapped system, the low doping behavior of fi is 
delicate. At T — 0, \i varies discontinuously at S = 0. 
This discontinuity stems from the energy gap between the 
valence and conductions bands. From both sides of the 
carrier doping close to 5 = 0, fj, is not a monotonic func- 
tion of 5. In contrast to this zero T behavior, at finite T, 
fx increases dramatically at small doping. It is well known 
that n at finite T is lower than that at T = for a one- 
band electron system with parabolic energy-momentum 
dispersion relation because the particles occupy energy 
levels above the Fermi energy with a lowered chemical 
potential. In the present gapped case, the conduction 
band is the effective one band system. At smaller elec- 
tron doping, the Fermi energy is lower, the temperature 
effect is much more pronounced. The situation for the 
zero-gap system is different. The temperature effect in /i 
of the zero-gap system is not so notable as in the gapped 
system. For A = 0, the electrons in the valence band 
can be thermally excited to high levels in the conduction 
band without significant change in the chemical potential 
because the density of states is nearly symmetric about 
the touch point of the valence and conduction bands. The 
chemical potential for A = at T/Tq = 1 is almost the 
same as at T = 0. 

The behavior of the chemical potential is closely re- 
lated to the compressibility k and the spin susceptibility 
X of the doped carrier system. They are denned as 



and 
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with [Ib as the Bohr magneton. The common factor 



[dnjd\x)T can be calculated by performing the deriva- 
tive with the obtained result as shown in Fig. 4. On the 
other hand, this factor is actually the irreducible density- 
density response function \nn diagrammatically shown in 
Fig. 2 (b) with the vertex correction given in Fig. 2 (c). 
Here, both results are the same because the SCHF ap- 
proximation for E(fc) and x„„ satisfies the microscopic 
conservation laws^l 

In Fig. 5, we show {dn/dn)i as a function of 5 for 
A = at various temperatures. At low doping close 
to 6 = 0, there is a sharp decrease in (dfx/dnjT at 
T = with decreasing S. But even close to S = 0, 
(dfi/dn)T is positive. This result is different from the 
existing calculation based on the perturbative HF ap- 
proach on a bilayer graphene models With increasing 
the temperature, {dfi/dn)T at zero doping increases and 
the function at low doping tends to show a flattened be- 
havior. Since (dfj,/dn)T > 0, the homogeneous electron 
system of A = is mechanically stable at any doping 
and any temperature. The result for the free electrons 
at T — is also plotted in Fig. 5, which is a monoton- 
ically decreasing function of \8\. Recent experiments on 
the compressibility^ seem qualitatively consistent with 
this free electron behavior. Why the Coulomb effect is 
not observed is still an open question. 

Shown in Fig. 6 is (dfj,/dn)T as a function of S for 
A/Eq = 0.02 at various temperatures. At T = 0, 
(dfi/dn)T is negative at small d consistent with the neg- 
ative slope of \x as shown in Fig. 4. This means that 
the homogeneous electron system with (dfi/dn)T < is 
unstable and favors to become an inhomogeneous state. 
At finite temperature, in principle, (dfj,/dn)T is positive 
at 5 = because of the behavior of fi as shown in Fig. 
4. At T > 0, it is seen from Fig. 6 that (d^/dri)? at 
S = is likely infinitive. At finite but low temperature, 
(dfi/dn)T decreases within a small region close to 6 = 0, 
then increases, and finally decreases slowly at large S. As 
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FIG. 6: (color online) {dfi/dn)T as function of doping concen- 
tration 8 for gapped system of A/Eo = 0.02 at T = (blue 
dash-dot line), T/To = 0.5 (green dashed line with squares), 
and T/To = 1 (red line with circles). 



shown in Fig. 6, at T/Tq = 1, (d[i/dn)T is positive at all 
the doping concentrations. We then deduce that the ho- 
mogeneous system is stable at any doping concentrations 
for T/T > 1. 

To know whether the system is stable or not is helpful 
for both the theoretical understanding of the thermody- 
namic properties and the technological application of bi- 
layer graphene. We have determined the phase diagram 
of the electron system with A/Eq = 0.02 as shown in 
Fig. 7. The red curve in Fig. 7 divides the T — 5 plane 
into homogeneous stable (k > 0) and unstable (k < 0) 
regions. Note that at the border, k and x g° to infinity. 
K — > oo means the system undergoes a phase separation 
or Wigner crystallization. % — >■ oo implies there is fer- 
romagnetization in the system. Therefore, in the region 
with (dn/dn)x < 0, the case of an inhomogeneous sys- 
tem, a Wigner crystal^ or a ferromagnetic stated may 
be possible. 

In summary, we have investigated the compressibility 
instability of the interacting electrons in bilayer graphene 



using the self-consistent Hartree-Fock approximation. 
We have studied the combined effects due to the Coulomb 
interactions and the energy gap between the valence and 
conduction bands in the chemical potential and the com- 
pressibility of the electrons. We find that the homoge- 
neous system with zero gap is always stable. For a gapped 
system, the compressibility becomes negative at low car- 
rier doping concentrations and low temperature, leading 
to the instability of the homogeneous system. The phase 
diagram distinguishing the stable and unstable regions 
of a typical gapped homogeneous system is given by the 
present calculation. 
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FIG. 7: (color online) T — 5 phase diagram of the typical 
gapped electron system of A/Eo = 0.02 in bilayer graphene. 
The red curve is the border between the homogeneous stable 
(ft > 0) and the unstable (ft < 0) states. 
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